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Abstract 

^^ ' A preferred form for the path integral discretization is suggested that allows the implementation 

Q^ , of canonical transformations in quantum theory. 

a^' 

^ ; 1 Canonical Transformations in Quantum Theory 

O, 

^L^ • Canonical transformations are of great utility in classical theory [Goldstein, 1980]; however they do not 

0^ ' work so well in quantum theory, giving rise to anomalous (order h ) potential like terms [Gervais and 

Jevicki, 1980; Klauder, 1980]. Despite the age of the problem of quantum canonical transformations, it 
^ . continues to inspire interest [Anderson, 1993; 1993; 1994; 1994; Swanson, 1994]. 

^y. • Although equally applicable to quantum field theory, there is no virtue in working within its added 

complexity, and the problem is discussed in the context of quantum mechanics. This issue is well in- 
vestigated in the setting of the Hamiltonian (phase space) path integral formulation of quantum theory 

0^ . [Feynman, 1948; 1965; Fanelli, 1976], due to its closeness to classical theory with its commuting variables. 



A transition amplitude is then generally written as: 



u 
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\^ . which is formal because it actually depends upon how it is discretized. In general, each integral of the 

5—( ' 

j^ ' discretized path integral has a leading error of order At, and this error is not lost in the limit of infinite 

time refinement because there are -^ integrals {T = ti, — ta). So this error makes a finite contribution in 

the final limit, and it is exactly in this sensitivity to At that the discretization scheme and equivalently 

the operator ordering of the operator formalism expresses itself [Schulman, 1981; Mayes and Dowker, 

1972; 1973]. However, if one could locate a scheme with error of higher order than At, no contribution 

would be hidden, and one might anticipate a better behaved object under formal manipulations. The 

idea is that in moving into this scheme from another, the hidden error would be exposed and correctly 

accommodated within the proposed canonical transformation. 

For an investigation of the error, end point integrations can be neglected, since the abandonment of 
any finite number of At corrections will make no difference to the total error in the small time limit. 
Advantage will be taken of this simplification throughout. 

Although there will be a formal investigation, it might be productive to first embark upon some 
numerical experiments to show the existence of the types of objects sought. 



2 A 'numerical' Investigation 

In the normal scheme, the standard path integral discretization is given by a variety of non-equivalent 
forms, a typical example with end points {qa,Pa), {<lb,Pb) being: 
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where At = ^ 
It is well known, and will be confirmed below (both numerically and analytically), that each of the 
above integrals has an error of order At. But just as the definition of the normal derivative has various 
forms of differing error (of relevance to computer numerics), namely: 

dx ,. x(t + At)-x(t) ,. 

—- = lim (3) 

dt At-o At ^ ' 

or 

dx ,. x(t ^ At) ~ x(t - At) 

—- = hm ^ —, — ^ (4) 

dt At^o At ^ ' 

the first being accurate to order At, and the second to order (At)^; so one might also seek a path 
integral scheme with higher order error. There is a candidate in the more symmetric (mid-point) form 
given by [Shiekh, 1988; Klauder, 1980; Daubechies and Klauder, 1985]: 



°? °? 2M 2M+1 , / . 2M , N \ 



where At — 



2M+1 

This symmetric discretization for the path integral exists only for even subdivisions, and the first 
three time subdivision refinements are listed in the appendix, where Wick rotation (t — > —it) has been 
performed to improve numerical convergence, having set h — 1 and dropped the end integrations, an act 
which does not effect the error contribution in the small time limit, but simplifies matters. 

The traditional form (equation 0) is known to have an error of order At, while the symmetric form 
might be anticipated to have a higher order error. This might first be investigated numerically before 
proceeding analytically, with all end points held to zero (a coherent state like path integral). 

These integrals differ from the limit by / = /^^ +ck(At)" -I- . . ., so (adopting the notation where the 
subscript on / indicates the number of q (or p) integrations): 

/ - / +«(^r)" + . . . (6) 
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Eliminating J^ and a in order to isolate the leading order error yields: 

/e - h 7-" - 5-" 



(9) 



U - k 5-" - 3-" 

For a leading order error of (At)^ one would get a ratio of | for small evolution times, while for (At)^ 
one would get a value of ||. Using the trial Hamiltonian of a simple harmonic oscillator, namely: 

H^\[p' + q^) (10) 

avoids the heavy numerical work involved in the accurate evaluation of a twelve dimensional integration 
(in fact, the integrals were performed analytically with a computer mathematics package). Using a short 
evolution time of T = 0.1, and fixing the end points at zero, leads to the result for the normal path 
integral: 

h-U .9978625573 ...- .9980047869 .. . _^_ ,.., 

— — .42829 . . . (11) 

U - h .9980047869 ... - .9983368743 ... ' ' 

the i (.42857...) confirming that the leading order error is indeed (At)^; while for the symmetric form 

one gets: 

Jfi-J4 .9951226146... -.9952203945... „^,,^ 

- .27557 . . . (12) 



.9952203945 ... - .9955752212 . 



the II (.27551...) corresponding to a leading order error of (At)^, so confirming the suspicion that a 
higher order scheme will be free of hidden contributions. 

3 Stochastic Terms 

The chance is taken here to derive the well-known results that in the path integral p ^ (At) ^2 and 
A<7 ^ (At)2 . An important exception to this rule is derived below. 

Beginning from the Hamiltonian path integral, which consists of many integrals of the form: 
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(13) 



It is assumed here that one is not dealing with unphysical, higher derivative theories (which would 
demand higher order 'momenta') so that the momentum is not found higher than quadratic order or 
negative powers. As a result the canonical transformation is somewhat limited in that it should not map 
to a theory of higher order, and the integral becomes: 



00 00 

J2 



JdqJ dpexp (1 (pAq- (^4^ +7(.,i)l> + V^(.,t)) At)) (14) 

— 00 —00 

The p integral may then be performed using the Gaussian resuh: 



e— '-/'^ds = J-e^ (15) 

V a 

Now since the p^ generates a 1/a, where a — iAt/2'mh, so each p (for even powers) contributes Hke 
(At)^2^ while p alone generates a /3^, i.e. — (7Ai — Aq)^/h^, so that p in odd powers contributes like 
(At)^. In performing the p integrals in equation n3 and so obtaining the Lagrange formalism; p becomes 
m{-^ —7)7 so that each Ag ~ (At) 2 [c.f. p ~ (At)^2] in even powers, and ^ (At)^ for odd powers. 
Another way to see this result is to expand out the p term, namely exp(ip(Ag — ^At)/h), and then note 
from symmetry that the contribution starts only at the second term, i.e. forp^(Aq — 7At)^, which indeed 
contributes like (At)^. This higher order contribution for odd powers will be crucial later. 

It is in this way that the contributing class of paths are seen to be stochastic (or Brownian) in nature. 
This behaviour of the path integral must be carefully accounted for when working to order At. 

4 Canonical Transformations in the Symmetric Path Integral 

In classical mechanics a canonical transformation is one that preserves the least action principle [Gold- 
stein, 1980]. For the path integral one might analogously require that there be a path integral represen- 
tation in the new variables (Q,P,t), if one existed in the old ones {q^p,t). The fact that the symmetric 
path integral has no 'hidden' parts should guarantee that formal canonical transformations are now valid. 
This in explicitly demonstrated below. 

A canonical transformation should be system independent, that is to say, the transformation should 
be canonical not only for some specific system, but for all problems with the same degrees of freedom. 
The amplitude may alter under such a transformation by at most a phase factor. So formally one gets: 

I Dq J -J^e^pl-Jipq~H{q,p,t))dt\ 

^00 -00 TT V *" / (16) 

= exp{j-{Ft-Fa)) J DQ J ^-^exp i-J i^PQ - K{Q,P,t)j dtj 

with F{q,Q,t) being an arbitrary smooth function. Since the above equation is to be true for all 
Hamiltonians, one gets: 

dF 
pq-H = PQ-K+— (17) 

the same condition as in classical mechanics, with F the generating function of the canonical trans- 
formation. For F — F{q, Q, t) one gets: 



dF 

pq-H = PQ-K+—- 

oq 



dF 
Q,t 9Q 



q.t 



dF 
Ik 



from which follows, by the independence of q and Q 

_ dF 
dq 



(18) 
q,Q 



(19) 



p= - 



K = H 



dF 
dQ 



q,t 



dF 

'at 



(20) 



(21) 
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All this work was formal, and now one is ready to apply this machinery to the symmetric discretization 
of the path integral in the hope that there will be no corrections. The formal canonical transformation 
of equation Hq now becomes (ignoring end integrals): 



oo 2M 
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A 2M 



Qk+1 — Qk-1 



^k=l' 



-K{Qk,Pk,tk)At + AF 



(22) 



where: AF = 



P{Qk+i, Qfc+i,ifc+i) — F{qk-i,Qk-i,tk-i) 



Now, F itself, like all other terms in the action, cannot be worse that (At)'^ in strength, for the same 
reason given before that 'higher derivative' actions are excluded from this discussion. By expanding F 
and also using the facts that Aq ^ (At) 2 and AQ ^ (At) 2 for even powers, and ^ (At)^ for odd powers, 
one gets the crucial result: 



AF^'-l 
dQ 



q,t 



AQ+ — 
oq 



Q,t 



dF 



At + 0(At)^ 



(23) 
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Had one not used the symmetric scheme, anomalous (order fi ) terms would have entered here with 
the even powers of Aq and AQ then presentjj. Dropping the (At) 2 term in the above (it is a At term 
that does not disappear in the limit) leads to: 



dF 

AF^-PAQ+pAq+ -- 

dt 



At 



(24) 



and from this one gets the anomaly free transformation of pAq, namely: 



pAq = PAQ + AF 



dF 

'dt 



At 



(25) 



q,Q 



The Jacobian for the measure is unity, and the Hamiltonian conversion is equally trivial, since they 
are both local, and leads to: 



H = K 



dF 

'dt 



(26) 



q,Q 



Putting this all together leads to the sought after equality, namely: 

^It was the lack of a viscous term in our numerical experiment that lead to a (At)^ error, as opposed to the more general 
(At) 2 error. 
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having restored the end integrals (at the cost of no error in the limit). 

This completes the demonstration that the symmetric form canonically transforms cleanly. 

5 Corollaries 

The utility of having a canonically invariant prescription for the path integral extends beyond just canon- 
ical transformations. It permits a consistent (although not unique) quantization of a classical system 
[Chernoff, 1981; Kapoor, 1984; Dirac, 1925; 1958]. Such consistency is part of the way to making the 
path integral a well-defined object [Daubcchies and Klauder, 1985]. The phase space path integral also 
has a sensitivity to the order in which the p and q integrations are performed. This problem and its 
solution is discussed elsewhere [Shickh, 1990]. 

It should be emphasised that the symmetric path integral is favoured (not compelled) over other 
prescriptions, in that it generates no anomalous additional terms during the canonical transformation. 
Note, however, that to get into and from this scheme, stochastic {h ) terms appear. The virtue of no 
terms occurring during the transformation is that the Hamiltonian behaves well under the transformation, 
and is then (for example) trivialised by a Hamilton- Jacobi transformation. 
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7 Appendix: Explicit Discretizations 

The first three (Wick rotated) symmetric path integral discretizations with h — 1 and end points held 
fixed (coherent state like path integral) are given by: 



(4^) 



dqidqi / dp2dp2 cxp 
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dqi..dq4 / dpi..dp4exp 



dqi..dqQ I dpi-Ap^cxp 



i(Pi('72 -qa))-H{qi,pi,ti)T/3 
+ ^{p2iqb - qi))-Hiq2,P2,t2)T/3 



( ^{pi{q2-qa))-H{qi,pi,h)T/b \ 
+ ^P2iq3 - qi))-H{q2,P2,t2)T/5 

+ §(P3(94-92))-ff(?3,P3,i3)T/5 

V +^{p4iqb-q3))-Hiqi,p4,U)T/5 J 

( ^{pi{q2-qa))-H{qi,pi,h)T/7 \ 
+ ^ip2{q3 - qi))-H{q2,P2,t2)T/7 
+ ^{P3{q4 - q2))-H{q3,p3,t3)T/7 
+ t(P4(<?5 - q3))-H{qi,pi,U)T/7 
+ ^{p5{q6 - q4))-H{q5,P5,t5)T/7 

V +UP6{qb - q5))-H{qe,Pe,te)T/7 J 



(28) 



(29) 



(30) 



These might be compared against the corresponding, more usual, discretizations given by: 



h = 



h = 



{2nY 



dqidqi / dp2dp2Cxp 
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qa) 

-qi) 

-q2) 



H{qa,Pl)T/3 

-H{qi,p2)T/3 
-H{q2,Pb)T/3 
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dqi..dqQ / dpi..dpQCxp 



( ipi{qi-qa) - H{qa,Pi)T/5 ^ 
+ip2{q2 - qi) - H{qi,p2)T/5 
+ip^{q^ - 92) - H{q2,P3)T/5 
+ip4iq4 - qs) - H{q3,p4)T/5 

V +ipb{qb - qi) - H{q4,pb)T/5 J 

( iPi{qi-qa) - H(qa,Pi)T/7 \ 
+ip2{q2 -qi) - H{qi,p2)T/7 
+-iP3{q3 - 92) - H{q2,P3)T/7 
+ipA{qA - qs) - H{q3,p4)T/7 
+ipb{q5 - qi) - H{qi,p5)T/7 
+ip6{q6 - qb) - H{q5,p6)T/7 

V +Wb{qb - qe) - H{qe,pb)T/7 J 
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